Phonon spectrum and dynamical stability of a quantum degenerate Bose-Fermi 

mixture 
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We calculate the phonon excitation spectrum in a zero-temperature boson-fermion mixture. We 
show how the sound velocity changes due to the boson-fermion interaction and we determine the 
dynamical stability regime of a homogeneous mixture. We identify a resonant phonon- exchange 
interaction between the fermions as the physical mechanism leading to the instability. 
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New developments in the field of ultracold atoms are 
creating opportunities for detailed studies of quantum 
degenerate Fermi gases [[[JD • The Pauli exclusion princi- 
ple results in the lack of s-wave scattering and the sup- 
pression of higher order scattering processes at low tem- 
peratures, so that the direct evaporative cooling spin- 
polarized fermionic atoms does not work. To circumvent 
this obstacle, one may either cool a Fermi gas containing 
two different spin states or mix the fermions with a gas 
of bosonic atoms and sympathetically cool this mixture 
||. Both methods have been successfully implemented 
experimentally to create degenerate Fermi gases. The ex- 
perimental progress along the latter direction has stim- 
ulated considerable renewed interest in studying Bose- 
Fermi mixtures in the quantum degenerate regime |3] ^) . 

The study of ultracold Bose-Fermi mixtures dates back 
to the 1960s, when the discovery that 3 He is miscible in 
4 He prompted the study of dilute solutions of 3 He in su- 
perfluid 4 He H. In that system, however, the physical 
picture is complicated by the strong interatomic inter- 
actions. As such, weakly-interacting alkali atomic vapor 
quantum gases provide us much cleaner systems to study 
the properties of Bose-Fermi mixtures. In this paper, we 
analyze the spectrum of bosonic elementary excitations 
in a Bose-Fermi mixture at zero temperature. In the 
case of an homogeneous system, we identify the stability 
regime of the mixture and calculate the sound velocity of 
the bosonic phonons under the influence of the fermions. 
We identify the physical mechanism leading to the in- 
stability of the mixture as a resonant phonon-exchange 
coupling between fermions below and above the Fermi 
sea. 

Our starting point is the Hamiltonian describing 
the interaction between a bosonic and spin-polarized 
fermionic system 
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'r^ (^ 0) -^ F )^F + f j #r$J> B fi F j> F , (1) 
where Ha^ = T a + V a (a = B, F) and the subscripts B 



and F represent 'boson' and 'fermion', respectively; T a , 
V a and n a are the associated kinetic energy, trapping po- 
tential, and chemical potential, and g and / represent the 
boson-boson and boson-fermion interaction strengths. In 
terms of the s-wave scattering lengths a, they are given 
by 

g = AttK cle/tub, f = I'kTi' o,bf lm r , 

with m r = rriBmp / (tub + mp) being the reduced mass. 
All higher order scattering processes are neglected. 

To determine the excitation spectrum of the Bose con- 
densate, we decompose the field operators as 



yl B {r,t) = 4> b {y) +£ fl (r,t), 
ip F (r, t)tp F (r, t) = n F (r) + 5p F (r, t), 



(2) 



where 0s(r) = (^(r, £)) is the condensate wave func- 
tion and n F (r) = (^(r, t)yli F (y, t)) is fermion density. 
Clearly, cf>B satisfies the Gross-Pitaevskii equation 



h b + 9 n B + fn F ) 4> B = 0, 



(3) 



where ns — |0s| 2 - 

Inserting Eqs. (|^) into (Q), neglecting terms involving 
three or more fluctuation operators yields for £b (r, t) and 
t) the Heisenberg equations of motion: 



ifi 



^B^B 



dt 



£l + I^bSpf, 
(H F - fi F )yj F + f(<f> B £,B + 4>*b^b)^f, 



(4) 
(5) 



with Hb = H^ 1 + 2griB + fnF and Hp = H 



(0) 



fn B - 



We first find an approximate solution of (j^). Ex- 
panding ijjp(r,t) on the basis of eigenstates <p n (r), with 
eigenenergies E n , of the Hamiltonian Hp as: 



$F(r,t) = a n (t)ip n (r), 



(6) 



and inserting this expansion into Eq. (||) yields the for- 
mal solution 
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xS(r',i> f (r',C), 



(7) 



where 



^ O) (r,f)=^a„(0)e-^ i /Vn(r) ) 



In first order, we replace tp F (r',t') by tjj F \r',t') on the 
r.h.s. of Eq. (]?]), so that 
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xS(r',f')f(r',0. 
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After some algebra, this expression yields the lowest or- 
der correction to the fermionic density in the presence of 
the Bose condensate as: 

Spf = i/J F (r,t)ip F (r,t) - n F (r) 

« 4 J dt> J dV J( - r ' r '' * ~ *') S ( r '' 

where 

,7(r, r', i - t') = G*> (r, r', t - t')G < (r, v',t- t') 

-G > (r,r',t-t')G* < (r,r',t-<'), 
G > (r,r',f-t')= E 



-iE n (t-t')/h * ( i 



^n( r 'Vn( r ): 
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and G< is the same as G>, but for £J„ < fj,p. In deriving 
Eq. (^), we have neglected second or higher order terms 
in £g and £ g , and we have replaced terms of the form 

$p "ftp by their expectation values with respect to the 
ground state. For example: 

i>P\r,t)ft°\v,t) « (i>^\r,t)^\v,t)) = n F (r). 

Inserting Eq. (|^) back into Eq. (|J) yields the equation of 
motion for the boson fluctuation operator: 

**^= (h { b } + 2gn B + fn F ) | B + g^ B 



results in a back-action onto the bosons — the last term 
at the r.h.s. of ©. 

Consider for concreteness a homogeneous mixture of 
Nb bosons and Np fermions confined in a box of vol- 
ume V. In the ground state, we have /ib = QTIb + f n F 
and hp — h 2 k F / (2m p) + fns , where we have introduced 
the densities ub = Nb/V and n F — Np/V = hp/ (<ott 2 ) 
and kp is the Fermi wave number. For periodic bound- 
ary conditions, the eigenstates of H F are plane waves 
tfk(r) = e tk ' r /\/V, with eigenenergies — h 2 k 2 /2m F + 
friB and k = |k|, so that the function J(y,t' ,t — t') of 
Eq. ( |l0| ) can now be explicitly expressed as: 



^(ry,t-o = ^ E E 

k n >kp k m <kf 



(A&-*£j(t-t') 



i(k m — k^Vr — i(k m — k n ) -r' 



Expanding likewise £s(r,i) as: 

CB(r,<) = -^E^We lk -' 



yields from Eq. (nO 
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at 



£ B (fc)r7k + ffn S 7?L k + (i/h)f 2 n B I(k) 7 (11) 



where Cs(k) = h k /2ms + gns, 

i"(k) = ~ /" tft' / dV /" d 3 re _ik,r J'(r,r',t-*') 

xE e ' ik '' r ' k*C) +%'(*')" 

k' 



c.c. 



k m <kf 



x6(|k m + k|-fc F ) »jl k (i') 



and O(x) is the unit step function. Similarly, 
lh ^T- = -£fl(*)»?Ik - gneVk - (i/h)f 2 n B I(k). (12) 



The two coupled integro-differential equations ( |ll|) and 
( |l2| ) can be solved using Laplace transformation. De- 
noting Ok (a) = I[r/k(t)] = J °° dte~ st T]k(t) and /3 k (s) = 

E( i s - (0) and © y icld 



+i ^L d *'y dVjr ( r ' r/ '*-* , ) 5 ( r/ '* / )^s(r). (10) 

This is one of the central results of this paper. It shows 
that the twofold effects of fermions on the excitation of 
the condensate: First, the fermions change the effective 
mean-field potential of the boson — the fn F term; second, 
the boson fluctuations affect the fermions, which then 



ihsa k (s) = C B (k)a k (s) + gn B Pk(s) 
f 2 

+i J - T n B ls(k){a^(s) + , (13a) 
n 

ifrs(3 k (s) = -C B (k)f3 k (s) - gn B au(s) 

-i^n B 4(A0M«) +&(*)], (13b) 



where 
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c.c. 



where ?iAk+k m ,k m = £ k m +k — £k m with ek = ft k /2mp. 
The excitation frequency uj is simply obtained by replac- 
ing s by iuj + + in Eqs. (|i~3|). Since these equations are 
homogeneous equations in ctk(s) and flk(s), uj is deter- 
mined by requiring the determinant of their coefficients 
to be equal to zero. These steps result in the phonon 
dispersion relation: 
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(14) 



e(|k n 



Ak+k m ,k„ 

k|-fc F ) 



uj + A k+km ,k m - i& 



-i0+ 
&{k F 



(15) 



Here we have replaced the sum over k m by an integral. 
As expected the Bogoliubov spectrum of a pure conden- 
sate is recovered for / = 0. 

Equation ([y|) has the same form as Eq. (12.30) in 
Ref. S. In general, i w {k) — and hence uj — may take com- 
plex values, indicating that the homogeneous state of the 
Bose-Fermi mixture becomes unstable. Its real part is 
given by the principal value of the integral in (Uq) and 
can be expressed as: 
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uj + a(k F + k/2) 



with a = ftk/m F . Its imaginary part, which determines 
the stability of the mixture, is given by 



Im[£ u (k)] =njdk m e(|k m + k| - kp)Q{k F - km) 

x <5(cj - A k+ k m ,k m ) , (16) 

which is nothing but the dynamical structure factor of 
the Fermi gas, S(k,uj), apart from a constant factor B. 
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FIG. 1. Resonant phonon-exchange interactions between 
fermions results in dynamical instability of the system. 



The physical interpretation of the instability associ- 
ated with Imfltpffc)] 7^ [or S(k,uj) ^ 0] is directly ap- 
parent from Eq. (|16[). The condition that a phonon of 
wave vector k and frequency uj satisfies u> = Ak+k m ,k m 
for some k m < k F and |k m + k| > kp implies that it res- 
onantly couples a pair of fermions lying inside and above 
the Fermi sea, respectively, see Fig. [y. When this cou- 
pling occurs, the Bose-Fermi mixture becomes dynami- 
cally unstable, and the system will spontaneously break 
into several distinct phases H. For long- wavelength 
phonons (k <ti kp), the stability criterion determined by 
Imjf^fc)] = requires the phonon frequency to satisfy: 



ftk 2 ftkkp 
uj > 1 . 



2mp 



mp 



(17) 
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FIG. 2. Stability diagram of a Bose-Fermi mixture. We 
assume mp = raj = m. We have adopted a system units 
in which the units for frequency, length, and wavenumber are 
hkp/2m, X/kp, and kp, respectively. Note that in this units, 
the fermion density is given by n F = 1/(6tt 2 ) « 0.017. 



The stability regime and the phonon spectrum can be 
obtained by first solving Eq. ( [Til) numerically while as- 
suming Im[^ w (fe)] = 0, and then checking if the system is 
stable using the criterion ([T7|). Fig. || shows the stability 
diagram of the mixture in the as F — nB space. As can be 
seen, the dynamical stability of the system is determined 
by both the scattering lengths and the densities. All 
other parameters being fixed, the stability condition (|l7 ) 
imposes a minimum boson density n™ beyond which 
no stable homogeneous mixture exists. We can use the 
Bogoliubov spectrum of a pure condensate to estimate 



nn mm 
9 n B 



ft 2 k 2 F 



v/rns 



ft 2 (6n 2 n F ) 2 / 3 
m F /m,B 



For realistic numbers, ji™" is about two orders of mag- 
nitude larger than n F (see Fig. ||). 
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FIG. 3. Phonon spectrum of a boson-fermion mixture. The 
dashed line corresponds to uj = h(k 2 +2kkF)/2rriF- The gray 
area below the dashed line represents the unstable regime. 
Same units as in Fig. 0. 



Figure || illustrates the phonon spectrum in the Bose- 
Fermi mixture for ub > n™. For small values of cibf, 
the spectrum is stable and single valued. Furthermore, 
the boson-fermion interaction increases the sound veloc- 
ity of the phonon. By expanding £ 0J (k) around small k, 
it is not difficult to show that in our dimensionless units 
(see caption of Fig. ^ and for raj = mp the sound ve- 
locity is approximately given by 



Co 



(2 + c ) 2 Tra B 



(18) 



where cq = \/16ira,BnB is the sound velocity in a pure 
condensate. For the parameters in Fig. we have a 
6% increase in sound velocity if clbf changes from to 
0.2. However, further increasing asF beyond a critical 
value Og F x splits the phonon spectrum into two branches, 
and one of them falls into the unstable regime. Again 
using the expansion of £ u (k), one can show that a B l p x is 
determined by (in dimensionless units): 



shown that a resonant phonon-exchange interaction be- 
tween fermions can induce a dynamical instability in the 
mixture. As an example, we have calculated the phonon 
spectrum of a homogeneous system and obtained analyti- 
cally the sound velocity under the effect of boson-fermion 
coupling. From the spectrum, we are able to determine 
the dynamical stability regime of such a system which de- 
pends on both the scattering lengths and the densities, 
and is directly linked to the dynamical structure factor 
of the Fermi gas. Our calculations show that for a given 
fcrmion density, a much denser boson field is required to 
stabilize the mixture. Viverit et al. also considered the 
stability of a homogeneous mixture In contrast to 
our work, which depicts a full dynamic picture for the 
boson-fermion mixture, their interest was focused on the 
mechanical stability phase diagram of the system, which 
is obtained by requiring the total energy to be a mini- 
mum but does not take the phonon-exchange interaction 
into full account. 

Our method can be straightforwardly applied to an 
inhomogeneous system, although the algebra in this case 
is much more involved. We will address this situation in 
a future publication. Another interesting case is when a 
boson field is mixed with a fermion field containing two 
components (e.g., two spin states). Interactions between 
the bosons and fermions will induce an effective coupling 
between the fermions from different components, which 
may result in the creation of Cooper pairs. 
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(clbf ) 2 — 27ra_BTO r /(ms + mp). 

This condition agrees with that obtained by Viverit et al. 
H and is confirmed by our numerical calculations (see 
Fig. |^) . For mp « raj , we have 



„max /ITI To 

a BF = V^a B /^- 



(19) 



Note that while n™" is independent of clef, aJ§ F x is 
independent of n B - Putting ( |l9|) into (|l8|), we obtain the 
maximum sound velocity achievable in a homogeneous 
mixture as: 



= c 
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In conclusion, we have developed a general formalism 
to determine the phonon spectrum of the condensate in 
a Bosc-Fcrmi mixture at zero temperature. We have 
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